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Abstract 

The glutathione transferase (GST) superfamily plays key roles in the detoxification of various xenobiotics. Here, we report 
the isolation and characterization of a silkworm protein belonging to a previously reported theta-class GST family. The 
enzyme (bmGSTT) catalyzes the reaction of glutathione with 1-chloro-2,4-dinitrobenzene, 1,2-epoxy-3-(4-nitrophenoxy)- 
propane, and 4-nitrophenethyl bromide. Mutagenesis of highly conserved residues in the catalytic site revealed that Glu66 
and Ser67 are important for enzymatic function. These results provide insights into the catalysis of glutathione conjugation 
in silkworm by bmGSTT and into the metabolism of exogenous chemical agents. 
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Introduction 

Glutathione (GSH) conjugation is essential for the detoxification 
of xenobiotics [1,2]. Several studies have also implicated 
conjugation reactions with endogenous compounds, such as a,(3- 
unsaturated aldehydes and prostaglandin [2-4], resulting in the 
excretion of at least one water-soluble compound. GSH transfer- 
ases (GSTs, EC 2.5.1.18) are responsible for catalysis of this 
conjugation and are distributed ubiquitously among aerobic 
organisms [5]. GSTs are cytosolic enzymes, widely distributed 
across both prokaryotic and eukaryotic kingdoms [6]. In 
mammals, there are seven GST classes (alpha, mu, pi, omega, 
sigma, theta, and zeta) that can be distinguished based on their 
primary amino acid sequence; identity is approximately 50% 
within a class and less that 30% between different classes [7,8]. Six 
GST classes (delta, epsilon, omega, sigma, theta, and zeta) have 
been identified in dipteran insects, such as Anopheles gambiae [9] and 
Dro.sophila melanogaster [10,11]. Insect GSTs can determine sensi- 
tivity to insecticides [9,12], and since the Lepidoptera are the 
principal insect pests in agriculture, knowledge of lepidopteran 
GSTs is of great importance. We have previously characterized 
several GSTs in the silkworm, Bombyx mori, a lepidopteran model 
insect [13—19], and a sigma-class GST in the fall webworm, 
Hyphantria cmea, one of the most serious lepidopteran pests of 
broad-leaved trees [16]. However, there have been no reports to 
date on the characterization of theta-class GSTs from silkworms. 

Here, we report the identification and classification of a theta- 
class GST isolated from B. mori, which we named bmGSTT. 
While bmGSTT shares some common substrates with human 
theta-class GSTs (hGSTT), it has a distinct substrate profile when 
compared to other B. mori GSTs studied to date. Furthermore, 
bmGSTT does not participate in the response to agents that 
generate oxidative stress, in contrast to previously identified B. mori 
GSTs. The activity profile of bmGSTT sheds further light on the 



way in which insects deal with xenobiotic agents and contributes to 
a more detailed understanding of the GST system in general. 

Materials and Methods 

Insects and tissue dissection 

Larvae of the silkworm, B. mori, were reared on mulberry leaves 
in the Institute of Genetic Resources, Kyushu University Graduate 
School (Fukuoka, Japan). At day -1 fifth instar larvae, fat bodies 
were dissected from the larvae on ice and stored at — 80°C until 
use. Total RNA was extracted rapidly from the dissected fat bodies 
with the RNeasy Plus Mini Kit (Qiagen Inc., Valencia, CA), in 
accordance with the manufacturer's instructions, and the resultant 
RNAs were subjected to RT-PCR. 

Cloning and sequencing of cDNA encoding bmGSTT 

Total RNA was processed using RT-PCR. First-strand cDNA 
was produced using Superscript II Reverse Transcriptase (Life 
Technologies, Carlsbad, CA) and an oligo-dT primer. The 
resulting cDNA was used as a PCR template with the oligonu- 
cleotide primers 5'-TATACCATGGTTTTAAAACTATATTA- 
TGAT-3' (sense) and 5'-CC GGATCC TTAAAGTTTAGAAT- 
TAGCCGCA-3' (antisense), based on a sequence obtained from 
the SilkBase EST database [20]. Underlined and double- 
underlined regions in the primer sequences represent jVcol and 
BamYLI restriction enzyme sites, respectively, which were used to 
insert the PCR product into an expression plasmid. PCR was 
performed with 1 cycle at 94°C for 2 min; then 35 cycles at 94°C 
for 1 min, 50°C for 1 min, and 72°C for 2 min; followed by 1 
cycle at 72°C for 10 min. The resulting bmGSTT cDNA (bmgstt) 
was ligated into the pGEM-T Easy Vector (Promega, Madison, 
WI), which was then used to transform E. colt DH5ot cells. Genetyx 
software (ver. 14.0.12, Genetyx Corp., Tokyo, Japan) was used to 
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bmGSTT 

human GSTT1-1 

human GSTT2-2 

AF515525 

AF515526 

CG1681 

CG1702 

CG30000 

CG30005 



pi al P2 a2a a2b |i3 

1 — mvlklyydlmsqpsrvlyillktmkydfepkyvnlrkaehysedftk-vnrmqrvpvi 5 7 

1 — mglelyldllsqpcravyifakkndipfelrivdlikgqIlsdafaq-vnplkkvpal 5 7 

1 — mglelfldlvsqpsravyifakkngiplelrtvdlvkgohIskeflq-inslgklptl 5 7 

1 msrsvklyydlmsqpsralyiflstnkipfdrcpialrkmqhkt: 'e : : rqvnrygkvpci 60 

1 msknlkyyydlmsqpsralwiflektklpyekclinlgkgehlteefk-ainrfqkvpci 5 9 

1 msqplkfyfdflnqssralyilleaskipfeaipismlkgehltgefrdnvnrfrklpai 60 

1 msapiryyydlmsqpsralfi ifrlsnmpfedcwalrngehltedfkkeinrfqrvpci 60 

1 mskaikyyydflsqpsralwiamklgktpfedcpvalrkqeqltdeyrs-inrfqkvpai 5 9 

1 mskpirfyydllspiarglwiglkfsnspveycpialrkfeqltdeykk-inrfqkvpai 59 

(54 a3 a4a a4b 

58 dhngfiltesiailkylsrenviaeslyskesklqarieeflewqhiglrlhcamyfrw 117 

58 kdgdftltBvaillyltrkykvpdywypqdlqararvdeylawqh|tlrrsclralwhk 117 

58 kdgdfiltBsailiylsckyqtpdhwypsdlqararvheylgwha|ci|gtfgiplwvq 117 

61 vdgsfrlaesvaiyrylcrefptdghwypsdtvrqarvdeylswqhlnlradvslyffhv 120 

60 tdsqiklaesvaifrylcreyqvpdhwypadsrrqalvdeylewqhhntratcai yfqyv 119 

61 tdhgyqlsenvaifrhlareklvpehwyprrhlgrsrideylawqqtnmgvatteyfqqk 120 
61 hdngyklaesvailrylsakgkipehlypkyfvdqsrvdeflewqhmslrltcamyfrtv 120 
60 vdgkfqlgesvsivryladkgvfseqlypktleerarvdeflewqhfnvrlvcslffrqv 119 
60 vggdfhlsetiai iryladkgqfdeklypktlenrarvdeflewqhlnirlacsmyfrda 119 



bmGSTT 

human GSTT1-1 

human GSTT2-2 

AF515525 

AF515526 

CG1681 

CG1702 

CG30000 

CG30005 



a4c a5 a6 

118 HMDPIL-TGRKSDEKTIQGYKRRMMMALDDFDTKWL — GRGTAFIVGETPTVADLVAACE 174 

118 VMFP-VFLGEPVSPQTLAATLAELDVTLQLLEDKFL QNKAFLTGPH I SLADLVAITE 173 

118 VLGP-L-IGVQVPEEKVERNRTAMDQALQWLEDKFL GDRPFLAGQQVTLADLMALEE 172 

121 WLNP-L-LGKEPDAGKTERLRRRLDGVLNFFDQELLSAGSGQAFLAGDRI S I ADLSAACE 178 

120 WLRPRM-FGTKVDPKQAEKYRGQMEGTLDFIEREYL — GSGARFIAGDE ITVADLLAACE 176 

121 WLVPYL-QKTRPADNAVNLASKQLEHTLNEFEQLFL NSRKFMMGDNISYADLSAICE 176 

121 WLEPLLT-GRTPSEAKIETFRMQMERNLDVVEEVWL EGKDFLTGSSLTVADIFAACE 176 

120 WLLPAKGLAPAPKPESVKKLIKDVESNLGLLERLWL E-KDFLVGDKLTVADIFGSSE 175 

12 0 WLFPMNGIAPKPKPEQIQALIEGVENNLGLLERLWL E-NDFLVGKNLTI-1ADILGSSE 175 



a7 a8 

bmGSTT 175 LEQPRMAGFEPKDH— FPNIAAWWPKVRDHFAPHYEDAHVILNKI INKMDRAANSKL 229 

human GSTT1-1 174 LMHPVGAGCQVFEGR-PKLATWRQRVEAAVGEDLFQEAHEVILKAKDFPPADPTIKQKLM 232 

human GSTT2-2 173 LMQPVALGYELFEGR-PRLAAWRGRVEAFLGAELCQEAHSIILSILEQAAKKTLPTPSPE 231 

AF515525 179 IEQAKIAGYDPCEGR— PALASWLTAVRERTNPYYDEAHKYVYRLSPDHIVTPVVAEDE — 235 

AF515526 177 IEQPRMAGYDPCEGR-PNLTQWMARVRESTNPYYDQAHKLVNKFAQDTASKAKL 229 

CG1681 177 IDQPKSIGYNAFQNR-NKLARWYETVREELGPHYKEVLGEFEAKLKGSGSGQQQGVAQAV 235 

CG1702 177 I EQTRMAD YDVRI K- YPK I RAWLKRVRQSCNPY YDVAHEFV YK I SGTGPQAKL 228 

CG30000 176 INQMKLCQYNVNEKQFPKVAKWMERVRDATNPYYDEAHSFVYKTSQQAVKAKN 228 

CG30005 176 INQLRLCQYRVDEKKFPKWKWLERVRVSANPYHDEGLTFIDRKSKQSTAAKL 228 

a9 

bmGSTT 230 229 

human GSTT1-1 233 PWVLAMIR 240 

human GSTT2-2 232 AYQAMLLRIARIP 244 

AF515525 236 235 

AF515526 230 229 

CG1681 236 KQ 237 

CG1702 229 228 

CG30000 229 228 

CG30005 229 228 



Figure 1. Alignment of amino acid sequences of theta-class GSTs. Sequences of GSTs were obtained in this study or from Swiss-Prot 
database (http://expasy.org/sprot/) and flybase (http://flybase.org): B. mori (present study), hGSTTM (No. P30711), hGSTT2-2 (No. PCG30), A. 
gambiae (AF515525 and AF515526), and D. melanogaster (CG1681, CG1702, CG30000, and CG30005). The G-site is shaded in blue. The oi-helices and 
p-strands are indicated in red and green and labeled with a and p. 
doi:1 0.1 371 /journal.pone.0097740.g001 

obtain the complete sequence of bmgstt and to deduce its 
corresponding amino acid sequence. Homology alignment (Fig. 1) 
was performed using ClustalW (ver. 1.83), with 10 and 0.2 as the 
gap creation penalty and gap extension, respectively. A phyloge- 
netic tree was generated using neighbor-joining plot software 
(http: / / www-igbmc.u-strasbg.fr/Bioinfo/ ClustulX/Top.html). 

Overexpression and purification of recombinant protein 

The bmgstt clone was digested with Ncol and BamHI and 
subcloned into the expression vector pET-llb, which was then 
used to transform competent E. coli Rosetta (DE3) pLysS cells 



(Novagen, EMD Biosciences, Inc., Darmstadt, Germany). Cells 
were then cultured at 37°C in Luria-Bertani media containing 
100 ug/mL ampicillin. After cell density reached an ODgoo of 0.7, 
isopropyl-l-thio-P-D-galactoside was added at a final concentra- 
tion of 1 mil to induce recombinant protein production. The 
culture was further incubated for 3 h, and cells were harvested by 
centrifugation. Bacteria were resuspended in 20 mM Tris-HCl 
buffer (pH 8.0) containing 0.5 M NaCl, 4 mg/mL lysozyme, and 
1 mM phenylmethanesulfonyl fluoride, and cells were subsequent- 
ly disrupted by sonication. Unless otherwise stated, all operations 
for purification described below were conducted at 4°C. The 
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supernatant containing the recombinant protein was clarified by 
centrifugation at 10,000 xg for 15 min and subjected to ammo- 
nium sulfate fractionation. The pellet obtained by ammonium 
sulfate fractionation was resuspended in 20 mM Tris-HCl buffer, 
pH 8.5. After dialysis against the same buffer, samples were 
subjected to anion-exchange chromatography on a DEAE- 
Sepharose column (GE Healthcare Bio-Sciences AB, Uppsala, 
Sweden) and eluted with a linear gradient of 0-0.3 M NaCl. The 
enzyme-containing fractions, assayed as described below, were 
pooled, concentrated using a centrifugal filter (Millipore Corp., 
Billerica, MA), and applied to a Superdex 200 column (GE 
Healthcare Bio-Sciences, Buckinghamshire, UK) equilibrated with 
the same buffer, but with the addition of 0.2 M NaCl. The purity 
of the pooled material was analyzed by SDS-PAGE using a 15% 
polyacrylamide slab gel containing 0.1% SDS, according to the 
method of Laemmli [21]. Protein bands were visualized by 
Coomassie Brilliant Blue R250 staining, and protein concentra- 
tions were measured using a Protein Assay Kit (Bio-Rad 
Laboratories, Inc., Hercules, CA), with bovine serum albumin as 
a standard. 

Molecular modeling 

A structural model of bmGSTT was constructed by SWISS- 
MODEL (http://swissmodel.expasy.org) [22] using the amino 
acid sequence. The model showed a GMQE (Global Model 
Quality Estimation) score of 0.69 [23]. The construction of the 
bmGSTT model was based on the structure of hGSTTl-1 (PDB 
code: 2C3T). The secondary structure assignments were produced 
with DSSP [24] . The Superpose program [25] revealed structural 
homology between bmGSTT and hGSTTl-1 with a root-mean- 
square deviation of 2.412 A/214 residues for all atoms. Figures 
were prepared using Coot [26] and PyMOL (http://pymol. 
sourceforge.net). 

Site-directed mutagenesis 

Amino acid-substituted mutants of bmGSTT were constructed 
using the Quick-Change Site-Directed Mutagenesis Kit (Strata- 
gene Corp., La Jolla, CA), according to the manufacturer's 
recommendations. An expression plasmid containing bmgstt was 
used as a template, and full-length mutated cDNAs were verified 
by DNA sequencing. 

Measurements of enzyme activity 

GST activity was spectrophotometrically measured using 1- 
chloro-2,4-dinitrobenzene (CDNB) and 5 mM GSH as standard 
substrates [27]. Enzymatic activity was expressed as mol CDNB 
conjugated with GSH per min per mg of protein. Alternatively, 
other substrates listed in Table 1 were used instead of CDNB 
[28,29]. Kinetic parameters (K m and k CM ) were assessed with a 
nonlinear least-squares data fit under assay conditions with 
different substrate concentrations in the presence of 5 mM 
GSH. The kinetic parameters toward GSH were measured in 
the presence of 1 mM CDNB. Thermostability of bmGSTT was 
determined by pre-incubation of enzyme solutions at various 
temperatures for 30 min before a residual activity assay. The pH 
stability of bmGSTT was assessed by pre-incubation of enzyme 
solutions at various pH values at 4°C for 24 h before a residual 
activity assay. Optimal pH for bmGSTT activity was determined 
using citrate-phosphate-borate buffer at various pH values with a 
fixed ionic strength of 0.25. 



Insecticide metabolism assay 

The ability of bmGSTT to metabolize each insecticide was 
determined by high performance liquid chromatography (HPLC) 
[14]. Reaction mixtures (500 uL) contained 120 uM PM, 
bmGSTT (12 ug), and 5 mM GSH in 50 mM Tris-HCl buffer 
atpH 8.0. Dehydrochlorinase activity for l,l,l-trichloro-2,2-bis(4- 
chlorophenyfjethane (DDT) was assayed by incubating the 
purified bmGSTT with 0.1 mM DDT and 5 mM GSH in 
20 mM Tris buffer (pH 8.0) at 30°C for 2 h. DDT and its 
metabolites were analyzed by HPLC, as described below, 
according to a previous report [14,30]. 

Reaction mixtures were extracted with three 500 uL portions of 
ethyl acetate for analysis by HPLC. After removing ethyl acetate, 
the amounts of each insecticide were determined by HPLC. An 
HPLC instrument (Prominence, Shimadzu Corp., Kyoto, Japan) 
was fitted with a 250x4.6 mm Mightysil RP-18 column (Kanto 
Chemical Co., Inc., Tokyo, Japan) with a flow rate of 1.0 mL/min 
at 40°C. The mobile phases employed were methanol (MeOH)/ 
acetonitrile/H 2 0 (70/15/15, v/v/v), MeOH/0.1% acetic acid 
(70/30, v/v), and MeOH/0.1% acetic acid (85/15, v/v) for 
detection of DDT, chlorfenapyr (CP), and permethrin (PM), 
respectively. The concentrations of each insecticide were deter- 
mined from the corresponding peak areas identified by its 
authentic sample. 

Results 

Sequence of cDNA encoding bmGSTT 

We isolated a cDNA, bmgstt, from silkworm fat bodies, the 
nucleotide sequence of which is deposited in GenBank under 
Accession No. AB848737. A BLAST search (http://blast.ncbi. 
nlm.nih.gov/Blast.cgi) using the Swiss-Prot database (http:// 
expasy.org/sprot/) showed that the sequence corresponds to 
theta-class GSTs. The sequence contains an open reading frame of 
690 base pairs, encoding 229 amino acid residues (Fig. 1), and the 
deduced amino acid sequence shows 34% and 31% identities to 
hGSTTl-1 and hGSTT2-2, respectively. In A. gambiae, there are 
two isoforms of the theta class (GenBank under accession 
numbers: AF515525 and AF515526), whereas, in D. melanogaster, 
4 isoforms of the theta class are known (flybase (http:/ /flybase.org) 
under accession numbers: CG1681, CG1702, CG30000, and 
CG30005). The amino acid sequence of bmGSTT reveals 
identities of 39%, 45%, 30%, 47%, 35%, and 35% to AF15525, 
AF15526, CG1681, CG1702, CG30000, and CG30005, respec- 
tively (Fig. 1). Secondary structural elements of bmGSTT were 
predicted using SWISS-MODEL (http://swissmodel.expasy.org/ 
workspace/). Additionally, the DSSP program [24] revealed that 
the bmGSTT monomer includes 9 a-helices and 4 f3-strands 
(Fig. 1). In hGSTTs, there are additional oc-helices, in comparison 
to bmGSTT, which would be present as a2a, a4c, and oc9. The 
locations for the other a-helices and f3-strands are conserved 
among the three GSTTs (Fig. 1). Construction of a phylogenetic 
tree clustered bmGSTT in the same clade as all current theta 
members (Fig. 2). The theoretical molecular mass and isoelectric 
point (26,913 Da and 8.91, respectively) of bmGSTT are similar to 
those of zeta- and delta-class GSTs from B. mori (Table 2). 

Putative GSH-binding site (G-site) 

The G-site identified in hGSTTl-1 include His40, Val54, 
Lys53, Glu66, Ser67, and Thrl04; whereas in hGSTT2-2, they 
are Lys41, Leu54, Glu66, Ser67, Aspl04, and Argl07 [31]. 
Superimposition of modeled bmGSTT on hGSTTl-1 indicates 
that equivalent B. mori residues include His40, Arg53, Val54, 
Glu55, Ser67, and Ilel04 (Fig. 3A). In this model, the distance is 
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Table 1. Substrate specificity of bmGSTT. 



Substrate 


Concentration (mM) 


Activity ([imol/min/mg) 


Wavelength [nm] 


Ae (mM H cm M | 


CDNB 


1.0 


0.03 


340 


9.6 


EPNP 


1.0 


2.57 


260 


0.5 


4NBC 


1.0 


NA 


310 


1.9 


A MDR 


1 .0 


3.56 


3 1 0 


1 .2 


4HNE 


0.1 


NA 


224 


13.8 


ECA 


1.0 


NA 


270 


5.0 


4NPA 


1.0 


NA 


400 


8.3 


H 2 0 2 


0.2 


NA 


340 


-6.2 


PM 


0.25 


NA 






DDT 


0.1 


NA 






CP 


0.25 


NA 







Activity was measured at pH 8 in the presence of 5 mM GSH. Data are expressed as means of three independent experiments. NA represents no activity. Wavelength 
and Ae represent maximum wavelength of the absorption and molecular coefficient, respectively. — : not applicable. 
doi:1 0.1 371 /journal.pone.0097740.t001 



large between the oxygen atom in the cysteinyl region of S- 
hexylglutathione (GTX), an inhibitor of GST, and the side-chain 
of Arg53 of bmGSTT (3.54 A), and there is a large distance 
between the amino group in the y-glutamyl region of GTX and 
the side-chain of He 104 (8.79 A) of bmGSTT. In hGSTT2-2, 
GSH interacts with amino acid residues (Lys41, Leu54, Glu66, 
Ser67, Aspl04, and Argl07) (Fig. 3B), which are superimposed 
upon Tyr41, Val54, Glu66, Ser67, Ilel04, and Argl07 in 
bmGSTT. In Fig. 3B, GSH is far away from Tyr41 (3.41 A 
between oxygen atom in the glycine portion of GSH and the 
hydroxy-group of Tyr41) and He 104 (8.39 A between the amino- 
group in the y-glutamyl region of GSH and the side-chain of 
He 104) in bmGSTT. Taken together, the structural data indicate 
that five residues (His40, Val54, Glu66, Ser67, and Argl07) in 
bmGSTT participate in the interaction with GSH. 

Characterization of bmGSTT 

Bacterially produced bmGSTT was purified to homogeneity, 
yielding a single band in SDS-PAGE with a molecular size of 
approximately 26,000 Da (Fig. 4). This size is close to the 
estimated size based on the deduced amino acid sequence. 

The enzymatic properties of purified bmGSTT were studied 
using CDNB and GSH as substrates. The enzyme was stable at 
temperatures <50°C and retained >75% of its original activity 
over a pH range of 5-1 1, similar to other B. mori GSTs. The pH 
optimum of bmGSTT was 8.0, which is identical to the optima of 
epsilon-, delta-, and sigma-class B. mori GSTs. 

Various substrates were then used to profile the activity of 
bmGSTT. Table 1 shows that bmGSTT was active toward 
CDNB, l,2-epoxy-3-(4-nitrophenoxy)-propane (EPNP), and 4- 
nitrophenethyl bromide (4NPB), but it was not active toward 4- 
nitrobenzyl chloride (4NBC), 4-hydroxynonenal (4HNE), ethacry- 
nic acid (ECA), or 4-nitrophenyl acetate (4NPA). The activity of 
bmGSTT toward CDNB is higher than that of GSTT-1 from rats 
(0.01 |0,mol/min/mg), mice (0.01 umol/min/mg), or humans 
(0.005 umol/min/mg) [32]. The activity of bmGSTT toward 
EPNP and 4NPB is lower than those of rat and mouse GSTTs but 
higher than that of hGSTTl-1 [32]. 4HNE is a preferred substrate 
for delta bmGST and omega bmGST (bmGSTO) (Table 2), 
whereas bmGSTT was relatively inactive against 4HNE. GSH 
peroxidase (GPx) activity was not observed in bmGSTT, unlike 
that in bmGSTO and epsilon bmGST. Subsequent HPLC 



analysis revealed that bmGSTT was unable to recognize DDT, 
CP, and PM as substrates (Table 1). 

Amino acid residues involved in catalytic function 

Based on the G-site of hGSTTl-1 and hGSTT2-2, we 
identified His40, Val54, Glu66, Ser67, and Argl07, as the 
candidate G-site of bmGSTT (Figs. 1 and 3). To determine 
whether these residues are important for catalytic activity, we 
performed site-directed mutagenesis. The resulting mutants were 
named H40A, V54A, E66A, S67A, and R107A and were purified 
from E. coli clones (Fig. 4). Each preparation of mutant enzyme 
was present as a single band in SDS-PAGE. Since the activity of 
bmGSTT toward EPNP and 4NBC did not fit the Michaelis- 
Menten equation, we determined kinetic parameters with CDNB, 
GSH, and 4NPB and compared these parameters with those of the 
wild-type (WT) enzyme (Table 3). With CDNB as the substrate, 
the enzyme's ^ was 1.5 mM, which was 3.8-, 3. 1-, 2.2-, and 0.96- 
fold the value for unclassified, delta-, omega-, and sigma-class 
GSTs, respectively [13,15,16,18]. The R m values for V54A, E66A, 
and S67A were higher than that of WT. The k cal values for V54A, 
E66A, and S67A were higher, while the values for other mutants 
were lower, compared to that of WT. The k cal IK m values from 
E66A and S67A were 41% and 28% of that of WT, respectively, 
and no large differences in k cat /K m values were observed for H40A, 
V54A, and R107A. With GSH as the substrate, the K m values for 
V54A and E66A were 3.1 and 3.6 times that of WT, whereas no 
K m could be calculated for the S67A mutant. The k cat /K m value of 
S67A was undetectable, whereas that for E66A decreased by 54%; 
no marked changes in k cal /K m values were observed for H40A, 
V54A, and R107A. With 4NPB as the substrate, the k cal /K m values 
for H40A and R107A were 22% and 40% of that of WT, 
respectively; a similar value was observed for V54A. For E66A and 
S67A, we were unable to detect the kcat/K m value with 4NPB. In 
summary, the most distinctive features of this mutagenesis are the 
decreased k ca i/K m values toward CDNB, GSH, and 4NPB for 
S67A, compared to those of WT. These results suggest that the 
interaction between GSH and Ser67 of bmGSTT is crucial for the 
activity. 
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Figure 2. Phylogenetic analysis of amino acid sequences of GSTs. The phylogenetic tree was produced by neighbor-joining plot software 
using GST sequences obtained from NCBI (http://www.ncbi.nlm.nih.gov/), Swiss-Prot (http://expasy.org/sprot/), and flybase (http://flybase.org/); each 
entry contains the species name, GST class, and accession number. Numbers represent branch length; Ag, A. gambiae; Md, Musca domestica, Dm, 0. 
melanogaster; Ms, Manduca sexta; He, Hyphantria cunea; Hs, Homo sapiens; bm, 6. mod; Ae, Aedes aegypti; SI, Spodoptera litura; Px, Papilio xuthus; and 
numbers attached to nodes indicate bootstrap values. The unclassified GST group does not include GST class. 
doi:1 0.1 371 /journal.pone.0097740.g002 



Discussion 

Although many GSTs have been identified in B. mori, the theta 
class remains poorly understood. This is a critical gap in our 
knowledge, because understanding the metabolic profile of theta- 



class GSTs may provide novel insecticide-targeting strategies. 
According to the silkworm genome sequence, there could be 23 
homologs of GSTs: delta-class (4 isoforms), epsilon-class (8 
isoforms), omega-class (4 isoforms), sigma-class (2 isoforms), 
theta-class (1 isoform), zeta-class (2 isoforms), and unclassified (2 
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isoforms) GSTs. In the A. gambiae genome, the GST classes include 
delta-class (12 isoforms), epsilon-class (8 isoforms), omega-class (1 
isoform), sigma-class (1 isoform), theta-class (2 isoforms), zeta-class 
(1 isoform), and unclassified (3 isoforms) GSTs, whereas, in D. 
melanogastet ; the classes include delta-class (1 1 isoforms), epsilon- 
class (14 isoforms), omega-class (5 isoforms), sigma-class (1 
isoform), theta-class (4 isoforms), zeta-class (2 isoforms), and no 
unclassified (0 isoform) GSTs. The silkworm genome contains a 
single gene encoding a theta-class GST. Previously, we reported 
identification of one theta-class GST of B. mori [13], which has 
been recently reassigned to the delta class. Thus, the focus of this 
study was on a silkworm GST in the theta class, which had not 
been thoroughly investigated, in terms of molecular and 
biochemical properties. 

GSTs catalyze a broad range of reactions, and each family 
member has its own discrete substrate specificity. This character- 
istic is also true for B. mori GSTs (Table 2). bmGSTT possesses 
GSH-conjugation activities toward EPNP and 4NPB, a property 
shared with mammalian theta-class GSTs. In contrast to 
hGSTTl-1, bmGSTT was not reactive with 4NBC and H 2 0 2 , 
suggesting that the catalytic properties of the bmGSTT enzyme 
are unique. bmGSTT did not recognize 4HNE, a cytosolic 
product of lipid peroxidation [33], or H 2 0 2 as substrates, 
indicating that the enzyme is unlikely to participate in the 
response to oxidative stress. Intriguingly, although bmGSTT 
shares some substrate preferences with mammalian GSTTs, it 
appears to have very different substrate specificity compared to 
other B. mori GSTs. Epsilon-class GSTs of mosquito could be 
involved in resistance to DDT and pyrethroid insecticides [34,35]. 
This resistance is particularly relevant given that HPLC analyses 
revealed that bmGSTT was unable to degrade the insecticides 
tested, in contrast to the results with other B. mori GSTs. 

The GST amino acid sequence is divided into two regions, the 
N- and C-terminal domains [5], The N-terminal domain includes 
the G-site, and the C-terminal domain has a hydrophobic 
substrate-binding site (H-site). The sequence diversity of the H- 
site dictates substrate selectivity [5] ; moreover, this diversity likely 
explains the varied substrate specificity of B. mori GSTs, because 
there is considerable divergence between their C-terminal regions 
(alignments not shown). Our mutagenesis results suggest that 
residues Glu66 and Ser67 in bmGSTT play important roles in its 
catalytic functions. Notably, while mutation of His40 in bmGSTT 
did not alter the kinetics of catalysis, the equivalent residue in 
delta- and epsilon-class GSTs is critical for GSH binding [14,36]. 
The mutation to Val54 had a minor effect on enzyme catalysis. 
This result was expected, because the mutation affected the main 
chain of the residue that interacts with GSH and not the side 
chain. We assume that His40 and Argl07 are not entirely crucial 
for binding of GSH and, instead, play co-operative roles with 
other residues in the G-site of bmGSTT. Similar observations 
were reported for an unclassified GST of B. mori (bmGSTu) [37], 
in which the equivalent residue (His53) of bmGSTu interacts with 
pre-bound GSH, but the mutation of the His to Ala did not affect 
catalytic activity. 

As mentioned above, the diversity of amino acids at the N- and 
C-terminal binding domains of GST is associated with substrate 
selectivity. hGSTTl-1 contains an H-site formed by Leu7, Leu35, 
Ile36, His40, Leulll, Trpll5, Metll9, Phel23, Hisl76, Leu231, 
Trp234, Val235, and Met238 [32]. We found that only 3 of these 
1 3 residues were conserved in the H-site of bmGSTT, which may 
explain the difference in substrate specificity between bmGSTT 
and hGSTTl-1. Additionally, a C-terminal helix in theta-class 
GSTs and residue 234 in the amino acid sequence of hGSTTl-1 
play important roles in substrate specificity and catalysis, 
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Figure 3. Superimposed structures of bmGSTT with hGSTTI-1 (A) and HGSTT2-2 (B) showing amino acid residues of the G-site. 

Carbon atoms for bmGSTT, hGSTTs, and GTX/GSH are green, cyan, and magenta, respectively, except for the regions of oxygen (red), nitrogen (blue), 
and sulfur (yellow). Symbols of amino acid residues for bmGSTT and hGSTTs are shown in green and cyan, respectively. 
doi:1 0.1 371 /journal.pone.0097740.g003 



respectively [31,32] [38]. There is no corresponding region, 
including the residue at position 234, in bmGSTT (Fig. 1), which 
may explain why it exhibits lower activity than rat, mouse, and 
human theta-class GSTs [39]. 

Recently, the electron-sharing network that contributes to the 
catalytic activity of GST has been described [40,41]. Based on an 
amino acid residue at position 64 that is functionally conserved in 



the GST classes [40], this network can be divided into type I and 
II classes. The type I electron-sharing network is exemplified by 
delta-, theta-, omega-, and tau-class GSTs, which contain an 
acidic amino acid residue at position 64, whereas the type II 
network GSTs (alpha, mu, and pi classes) have a polar amino acid 
residue. Glu66 is conserved in the sequence of bmGSTT; thus, this 
enzyme resembles a member of the type I network. The electron- 
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Figure 4. Recombinant bmGSTT and its mutants overexpressed in £ coli. bmGSTT mutants were prepared by site-directed mutagenesis. 
Purified recombinant proteins were applied to SDS-PAGE followed by Coomassie Brilliant Blue R250 staining; M, protein molecular size markers; lane 
1, WT; lane 2, H40A; lane 3, E66A; lane 4, S67A; lane 5, R107A; and lane 6, V54A. 
doi:1 0.1 371 /journal.pone.0097740.g004 
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Table 3. Comparison of kinetic data from bmGSTT and mutant forms. 





bmGSTT mutants 




WT 


H40A 


V54A 


E66A 


S67A 


R107A 


CDNB 


Km 


1 .5 (0.25) 


0.52 (0.080) 


28 (3.6) 


12 (2.1) 


4.3 (0.72) 


0.96 (0.11) 


k b 


0.27 (0.020) 


0.12 (0.036) 


3.9 (0.84) 


0.88 (0.85) 


0.35 (0.022) 


0.16 (0.042) 




0.18 


0.23 


0.22 


0.073 


0.051 


0.17 


GSH 




3.6 (0.56) 


1.8 (0.34) 


11 (1.9) 


12.9 (2.2) 


ND 


3.9 (0.35) 


k b 


0.15 (0.034) 


0.065 (0.018) 


0.52 (0.082) 


0.28 (0.077) 


ND 


0.13 (0.022) 




0.041 


0.036 


0.049 


0.022 


ND 


0.033 


4NPB 


Km 


0.13 (0.03) 


2.3 (0.22) 


0.50 (0.083) 


ND 


ND 


0.52 (0.071) 


k b 


0.78 (0.17) 


3.1 (0.56) 


3.7 (0.61) 


ND 


ND 


1,3 (0.16) 




6.0 


1.3 


7.4 


ND 


ND 


2.4 



Values, except those of k cat /K m , are expressed as means of three independent experiments. 

Expressed in units of mM; b expressed in units of umol/L/min; and Expressed in units of/min/mM. ND represents 'not detected'. 
doi:1 0.1 371 /joumal.pone.0097740.t003 



sharing network in hGSTT2-2 was proposed to contain Ser67 as 
one of residues involved in the network [41]. The equivalent 
residue in bmGSTT (Ser67) is conserved (Fig. 1). Glu66 and Ser67 
in bmGSTT could be part of an electron-sharing network and the 
G-site via direct interaction with GSH. Thus, mutation of the 
residues may result in a decrease in GSH-conjugation activity. 

Other than five residues (His40, Val54, Glu66, Ser67, and 
Argl07) in bmGSTT, there could be other amino acid residues 
that are essential for bmGSTT catalytic activity. In theta-class 
GSTs, the Ser residue in the N-terminal domain is conserved [42— 
44] and considered important for activation of the bound GSH. 
The equivalent residue in bmGSTT is Serll (Fig. 1). In other 
GST classes, mutagenesis of amino acid residues in electron- 
sharing networks results in decreased activity [40,41], Investiga- 



tion of putative catalytic residues using site-directed mutagenesis is 
now underway in our laboratories. 

Our results suggest that bmGSTT might play a role in 
detoxification of xenobiotics in B. mori. Together with bmGSTT, 
the roles of other GSTs in B. mori should be further examined to 
understand the mechanisms underlying insecticide detoxification. 
In turn, such studies will aid the design and implementation of 
insecticide-resistance management strategies for agricultural pests. 
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